Thermoelectric phenomena in a quantum dot asymmetrically coupled to external leads 
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We study thermoelectric phenomena in a system consisting of strongly correlated quantum dot 
coupled to external leads in the Kondo regime. We calculate linear and nonlinear electrical and 
thermal conductance and thermopower of the quantum dot and discuss the role of asymmetry in 
the couplings to external electrodes. In the linear regime electrical and thermal conductances are 
modified, while thermopower remains unchanged. In the nonlinear regime the Kondo resonance 
in differential conductance develops at non-zero source-drain voltage, which has important conse- 
quences on thermoelectric properties of the system and the thermopower starts to depend on the 
asymmetry. We also discuss Wiedemann-Franz relation, thermoelectric figure of merit and validity 
of the Mott formula for thermopower. 
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PACS numbers: 75.20.Hr, 72.15.Qm, 72.25.-b 

I. INTRODUCTION 

The Kondo effect, discovered more than seventy years 
ago in metals with small amounts of magnetic impuri- 
ties and explained more than forty years ago [l|, is now 
a prime example of many body phenomena in a wide 
class of correlated electron systems [2j. Due to recent 
advances in nanotechnology it became possible to fab- 
ricate artificial structures, where the Kondo effect can 
be studied systematically. Probably the best control 
over the isolated spin can be achieved in quantum dot 
(QD) systems. Originally, the Kondo effect in quantum 
dots was predicted theoretically in late eighties [1, d, Q 
and later confirmed in series of beautiful experiments 

UMMMMMMMM- 

The origin of the Kondo effect, both in alloys and in 
quantum dots comes from the formation of a singlet state 
between a localized spin and free electrons. This takes 
place at temperatures lower than the Kondo temperature 
Tk and manifests itself in increasing conductance (G) of 
the system. The increase of G is a result of the formation 
of a many body Kondo or Abrikosov-Suhl resonance at 
the Fermi energy. 

The experiments have confirmed the validity of the 
theoretical picture but they also discovered phenomena 
which require new theoretical ideas. These include obser- 
vation of the Kondo resonance at non-zero source-drain 
voltage d, Q , absence of even-odd parity effects expected 
for these systems [T(|, observation of the singlet-triplet 
transition in a magnetic field 
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splitting of the Kondo 
resonance due to the ferromagnetism in the leads [l2j 
or the interplay between the Kondo effect and supercon- 
ductivity in carbon nanotube quantum dots coupled to 
superconducting leads [l3j . 

In the present paper we shall focus our attention on 
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the nonequilibrium Kondo effect (i.e. the Kondo effect 
at non-zero source drain voltage) in asymmetrically cou- 
pled quantum dots. In our previous work [l5| we have 
found that the asymmetry in the couplings to the leads 
is the main reason responsible for the nonequilibrium 
Kondo effect. The effects of non-symmetric coupling 
have also been discussed in Refs. [3, [I?], 0, [H, H3|- 
Here we present a more systematic study on the role 
of the asymmetry in the couplings to the leads in ther- 
moelectric and transport properties of quantum dot sys- 
tem. Thermoelectric properties (thermopower and ther- 
mal conductance) in the Kondo regime have already been 
investigated in the quantum dots coupled to the normal 
[2ll . l22l [23j and the ferromagnetic [U, [25| leads, how- 
ever the role of the asymmetry in the couplings has not 
been discussed so far. Thus we shall concentrate on the 
conductance, thermal conductance, thermopower and re- 
lated quantities such as thermoelectric figure of merit 
which directly provides the information on the usefulness 
of the system for applications and Wiedemann-Franz ra- 
tio which signals breakdown of the Fermi liquid state 
when its normalized value differs from 1. 



In order to calculate those quantities we use non- 
crossing approximation (NCA) [26|, which is a widely 
accepted and reliable technique to study the equilib- 
rium and nonequilibrium Kondo effect in quantum dots 
[lH I27I I2H [2^ | . Although this method is known to give 
non- Fermi liquid ground state, it remains valid down to 
temperatures below Tk 0, El • 



The paper is organized as follows. In Sec. [IT] we in- 
troduce the model and discuss some aspects of our pro- 
cedure. The results of our calculations regarding the 
electrical and thermal transport in linear and nonlinear 
regime are presented and discussed in Sees. [TTTJ and IIV1 
respectively. Finally, conclusions are given in Sec. [V] 
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II. THE MODEL AND APPROACH 

The system we are studying consists of a quantum dot 
with a single energy level coupled to external electrodes. 
We describe it by the single impurity Anderson model 
[32| with very strong on-dot Coulomb repulsion (U — > 00) 
and use slave boson representation [33j , in which the real 
on-dot electron operator d a is replaced by the product of 
boson b and fermion operators f a (d a = b + f a ) subject 
to the constraint b + b + fa f<? = 1- The resulting 
Hamiltonian reads 

h = E £Akc^ k(T c Ak(T + e d ^2 Ufa + 

E(^4 kCT 6 + / CT + H.c), (1) 

Ak 

where A = L (R) denotes left (right) lead, c^ ko . (c.\ker) 
is the creation (annihilation) operator for a conduction 
electron with the wave vector k, spin a in the lead A 
and V\k is the hybridization parameter between local- 
ized electron on the dot with energy Ed and conduction 
electron of energy e A k in the lead A. 

In order to calculate electrical current J e \ and energy 
flux Jex flowing from the lead A to the central region 
we follow standard derivation (34[ and express all the 
currents in terms of Keldysh Green functions [35[ . More- 
over, we use relation Jq\ = Jex — ^xJeX for thermal flux 
JqXj so the resulting expressions are 

/OO 7 
^ —T x (u)[G<(u) + 2 l / A ( W )Im^( W )](2) 

J ^ = ^E jT ^r A H( w -MA)[G<(«) 

+2i/ A (w)ImGS(w)], (3) 

where = 2ir ^ k |VAk| 2 <5(<^ — £Ak) is the cou- 

pling of the dot to the lead A, and G r a (u>) is the 
time Fourier transform of retarded Green function 
(GF) G£(t,f) = ie(t-t'){[b+(t)U(t), f+(t')b(t')}+) and 
G<(w) = i{ft{t')b{t')b + {t)fJt)) is the Fourier trans- 
form of lesser Keldysh GF j3j|. f\(u) = (exp(f^-) + 
1) _1 is the Fermi distribution function in the lead A with 
chemical potential [i\ and temperature T\. 

In order to calculate lesser GF G<(u) we use widely 
accepted Ng ansatz 36] , as in the presence of both the 
on-dot Coulomb interaction and tunneling between the 
QD and the leads, it is not possible to calculate G<(u>) 
exactly. In this approach one assumes that full interact- 
ing lesser self-energy is proportional to the noninteract- 
ing one, and the resulting lesser self-energy is expressed 
in terms of retarded interacting and noninteracting self- 
energies. This approach has three advantages, (i) it is 
exact in nonequilibrium for {7 = 0, (ii) it is exact in 
equilibrium for any U, and (iii) it satisfies the continu- 
ity equation Jl — —Jr in the steady state (36| . In the 



present case it yields 

= -f E C ^ f M!/iH - /«H]imG;H (4) 
^iEOrH(^) 

x[f L (cj) - f R (cj)}lmG r a (cj), (5) 

where T = T L T R /(T L + Tr), and eV — \x h - \x R . The 
on-dot retarded GF G r a (w) is calculated within NCA in 
terms of boson and fermion propagators, which are ex- 
pressed by the coupled integral equations [HI, H3, [H, Hi| . 
In the numerical calculations we have used Lorentzian 
bands in the electrodes of width D = 100r, and chosen 
r = + T R 1 = 1 as an energy unit. The asymmetry in 
the couplings to the leads is defined as 5 = Tl/Tr. Note 
that in our previous paper [l5l] , to faithfully reflect the 
situation in experiment we kept Tr = 1 and varied 
r^, so T^+Tr ^ l. Here we shall concentrate on scaling 
of various physical quantities with asymmetry parameter 
S, therefore we keep T = 1, which means that both cou- 
pling parameters and Tr change at the same time. 

In the linear regime, i.e. for small voltage biases 
(eV = fj-L ~ Mil — * 0) an d small temperature gradi- 
ents (ST = Tl — Tr — > 0) one defines the conductance 
G = -(e 2 /T)L n , thermopower S = -(l/eT)(L 12 /L n ) 
and thermal conductance k = (1/T 2 )(L22 — L\ 2 jL\\). 
The kinetic coefficients read 

^ = ?E/^Hi mG ;M(-^M) T (6) 




with equilibrium Fermi distribution function /l(^) = 
/«(«) = /(«). 

In the nonlinear regime, i.e. for any voltages and 
temperature differences, it is not possible to use ki- 
netic coefficients, as they have been derived for small 
deviations from equilibrium [37| . Therefore, in gen- 
eral case, the thermoelectric quantities have to be cal- 
culated from proper general definitions, like the nonlin- 
ear differential conductance G(eV) = dJ e /d(eV), ther- 
mopower S = AV/ AT\j c= q or thermal conductance 
K = -Jq/AT\ Js=0 . 

III. LINEAR REGIME 

In the linear regime the transport coefficients depend 
on the asymmetry parameter via T only. This stems from 
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the fact that the equilibrium density of states (DOS) en- 
tering Eqs. ([H])-© is independent of asymmetry as long 
as Tl + Tr = 1. As a result Lij(8) can be easily ob- 
tained from its behavior in the case of symmetric cou- 
plings (<5 = 1) as 



L ij {8)=AL\ 



(l + 8f 



(9) 



where by L*j we denoted the corresponding kinetic coef- 
ficient for symmetrically coupled QD, i.e. Lf = Li AS = 
1). 

Knowing 8 dependence of the Ly, one can deduce 
what will be the modifications of the transport prop- 
erties due to the asymmetry in the couplings. For ex- 
ample, the linear conductance is directly related to the 
kinetic coefficient L\\ via the previously mentioned rela- 
tion G = —(e 2 /T)Li 1 . Thus G(S) shows the same scaling 
with 8 as Lij docs. It is linear in 8 in the limit 8 << 1 
and goes as 1/8 for 8 » 1. 

We shall not present the numerical results for conduc- 
tance and thermal conductance as for 8 = 1 as they are 
well known from previous studies [111, H2, H3] • The same 
is true for thermopower. However in Fig. [1] we show 
the temperature dependence of the (8 independent) ther- 
mopower calculated from kinetic coefficients L^ (solid 
line) and compare it to that obtained from so called Mott 
formula [38| (dashed line). 




At this point we would like to comment on the validity 
of so-called Mott formula for thermopower [38| , which is 
widely used in the literature EH, IHJ , while explaining 
experimental results. This formula relates thermopower 
and the linear conductance, i.e. 
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7r 2 fc| T dG 
3e GdE^ 



(10) 



and is expected to work well in noninteracting systems 
at low temperatures [43[ . As we can see in Fig. Q] the 
Mott formula works surprisingly well in the interacting 
systems in the whole temperature range. The qualitative 
behavior of Sm is approximately the same as of linear 
thermopower 5* (solid line) , giving slightly different tem- 
perature at which Sm changes sign. It also does not de- 
pend on the asymmetry 8. Thus one can conclude that 
the interactions do not necessarily lead to the violation of 
the Mott formula. Similar conclusions have been recently 
obtained for interacting quantum wires (44j . 

The low temperature thermopower has been measured 
in the Kondo regime in quantum dot system by Scheib- 
ner et al. [4C| . These authors found departures from the 
Mott formula and attributed them to developement of 
the narrow Kondo resonance. On the other hand, the 
measurements of thermopower in carbon nanotubes [4lj 
and in Zn nanowires [421 ] qualitatively agree with the 
Mott formula, Eq. ([10]), 

Figure [5] shows temperature dependence of the 
Wiedemann-Franz (WF) law which relates thermal and 
electrical transport via relation 
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This law describes transport in Fermi liquid bulk metals 
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FIG. 1: The linear thermopower S (solid line) and the Mott 
thermopower Sm (dashed line) as a function of temperature. 
Both quantities are not sensitive to the asymmetry S. 
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Thermopower is very useful measure of the Kondo cor- 
relations [H], El- 0, IS] . Below the Kondo temperature 
Tk it is negative, indicating particle like transport and 
changes sign at higher temperatures, where the transport 
is hole-like. At T « V it shows a broad maximum associ- 
ated with single particle excitations. The sign change can 
also be understood from the fact that S is sensitive to the 
slope of the density of states at the Fermi energy. With 
decreasing temperature, the Kondo correlations lead to 
development of a narrow peak in the DOS slightly above 
the Fermi energy, and thus to a slope change at the Ep. 
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FIG. 2: The Wiedemann-Franz ratio 3e 2 /tt 2 (k/TG) as a func- 
tion of temperature. This quantity is not sensitive to the 
asymmetry 8. 

and in general is not obeyed in QD systems where the 
transport takes place through a small confined region 21. 
[22l . [2J|. However, at very low temperatures, where the 
Kondo effect develops and the ground state of the system 
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has Fermi liquid nature, the WF law is recovered. At high 
temperatures, the WF law is violated as the transport is 
due to sequential processes leading to the suppression 
of the thermal transport (24|. The WF relation docs 
not depend on the asymmetry 8, as both electrical and 
thermal conductances show similar scaling with 5. 

At this point we would like to comment on the valid- 
ity of NCA, as it is well known, that this approach gives 
non- Fermi liquid ground state |30l. l3l|. However, around 
the Kondo temperature and slightly below Tk it provides 
a proper description. In our case Tk = 5.3 • 10 _2 r, thus 
NCA is believed to give reliable results in the tempera- 
ture range studied here. However, at lower temperatures 
we expect a violation of WF law. 

Thermoelectric figure of merit Z — S 2 G/k is a direct 
measure of the usefulness of the system for applications. 
For simple systems it is inversely proportional to opera- 
tion temperature, thus it is more convenient to plot ZT, 
which numerical value is an indicator of the system per- 
formance. Figure [3] shows ZT as a function of temper- 
ature. It has two maxima, one associated with charge 
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fiR. To be precise, there will be two resonances, one 
located at u> = fiL and the other one located at u = /sr. 
If we introduce asymmetry in the couplings S, one of 
the resonances will be broader and higher and the other 
one will be narrower and lower, depending which of the 
couplings Tl or Tr is stronger [15j | . As a results, all 
the thermoelectric quantities will be modified in more 
complicated way. 

The most pronounced example is the differential con- 
ductance Gdiff = dJ e /d(eV). In symmetrically coupled 
QD there is a resonance at zero source-drain bias. The 
introduction of asymmetry in the couplings (8 ^= I) shifts 
the resonance to non-zero voltage [15| . Figure [4] shows 
differential conductance Gdiff for different values of the 
asymmetry parameter 6. The spectra are calculated for 
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FIG. 4: Differential conductance Gdiff as a function of the 
source-drain voltage eV for different coupling asymmetries 
8 = 1, 2, 5 and 10 (from top to bottom), calculated at tem- 
perature T — 10~~ 2 I\ Note evolution of the position of the 
resonance with 8. 



FIG. 3: Temperature dependence of the thermoelectric figure 
of merit ZT = S 2 GT/k. This quantity does not depend on 
the asymmetry 8. 

fluctuations around T = T and the other one associated 
with the Kondo correlations below the Kondo tempera- 
ture. At temperature at which thermopower S changes 
sign it vanishes. However, as one can see, its value never 
exceeds 1, which indicates limited practical applicability 
of the system for cryogenic purposes. Similarly as the 
Wiedemann-Franz relation, it also does not depend on 
the asymmetry in the couplings. 

IV. NONLINEAR TRANSPORT 

In nonlinear regime we do not expect simple scaling of 
thermoelectric quantities with the asymmetry 8, as the 8 
will modify not only the effective coupling T in Eqs. (|4]) 
and {5J but also the QD density of states. In general, 
under nonequilibrium conditions the Kondo resonance in 
the DOS will be split by the bias voltage eV — — 



temperature T = 10 _2 r, below the Kondo temperature 
Tk, which is equal to 5.3 • 10 _2 r in this case. One ob- 
serves the evolution of the position of the maximum to- 
wards negative voltages with increasing of the 8. This is 
the effect of different heights and widths of two Kondo 
resonances in the density of states, one at [ih and the 
other one at \xr [15]. It is also clearly seen that the asym- 
metry suppresses the conductance for all bias voltages. 
This is due to suppression of the effective coupling T in 
Eq. by the asymmetry parameter 8. At lower tem- 
peratures the maximum of the conductance is narrower 
due to temperature effects and similar suppression of the 
conductance with increasing 8 is observed. In Fig. [5l 
the position of the conductance maximum eVo is plotted 
as a function of the asymmetry parameter 8 for various 
temperatures. The reason for larger shifts of the max- 
ima at higher temperatures stems from two facts. First 
one is the asymmetry in the couplings, which leads to 
a broader and higher (narrower and lower) Kondo res- 
onance associated with tunneling between QD and left 
(right) lead. The second reason, more important, is the 
temperature effect. At very low temperature, the energy 
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FIG. 5: The position of maximum in differential conductance 
as a function of the asymmetry parameter 8 for various tem- 
peratures T. Each position of the maximum has been ex- 
tracted from corresponding dJ e /d(eV) vs. eV characteristics 
with uncertainty given by errorbars. 



window accessible to the transport is bounded by the 
positions of the chemical potentials fiL and So in 
fact, only half of each Kondo resonance contributes to 
transport. At higher temperatures, this energy window 
becomes wider, thus it increases effect of the asymmetry, 
leading to larger shifts of the maxima in Gdiff- All curves 
in Fig. [5]show f(S) = a + b/S behavior with temperature 
dependent coefficients a and b. 

Similarly, the maximal value of the differential conduc- 
tance (associated with the Kondo effect), i.e Gdiff(eVb), 
also depends on the asymmetry parameter 5, as shown 
in Fig. [6] 
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FIG. 6: Peak height in Gam vs - asymmetry parameter 8 for 
a number of temperatures T. 

Another quantity affected by the asymmetry in the 
couplings is the nonlinear thermopower, shown in Fig. [7] 
as a function of the left lead temperature Tl- The right 
lead is kept at fixed temperature equal to 10 _3 r. Simi- 
larly as in the linear regime, the thermopower is negative 
in the low temperature (Kondo) regime and positive at 
high temperatures. The asymmetry 8 leads to shallower 
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FIG. 7: Nonlinear thermopower S as a function of the left lead 
temperature Tl for various values of the asymmetry parame- 
ter 8. The temperature of the right electrode Tr is fixed and 
equal to 10 _3 r. Note that, unlike linear thermopower (Fig. 
[I]), the nonlinear S strongly depends on 8. Inset: Normalized 
thermopower (S(8)/S(S = 1) as a function of the asymmetry 
8 calculated for T L = 1.9 ■ 10~ 2 r. 



minimum of S in the Kondo regime, and to higher val- 
ues of S at temperature around T, associated with single 
particle excitations. The change of the minimum of S at 
low temperatures reflects a fact that asymmetry changes 
the slope of the QD density of states at the Fermi energy 
due to splitting of the Kondo resonance under nonequi- 
librium. There is no simple scaling behavior of nonlinear 
thermopower with asymmetry 5. 

Another interesting fact is that S changes sign at dif- 
ferent temperatures for different asymmetries 5, leading 
to different Kondo temperatures. At first sight this is 
in contradiction with the conclusions reached from the 
linear behavior of the thermopower. However one has 
to keep in mind that for the system out of equilibrium 
(Pl 7^ fJ-R and T L ^ Tr) one of the Kondo peaks may 
disappear. For the parameters studied we always have 
the Kondo effect associated with the right lead but not 
with the left one. The density of states has a Kondo res- 
onance at the energy lu — /ir, while the Kondo resonance 
at fi L vanishes with increasing T L . If the asymmetry in 
the couplings is introduced, the Kondo resonance at ^lr 
contributes less and less to electrical and thermal trans- 
port. As a result, for large asymmetry 5 one observes 
lower minimum at low Tl and larger values of S at higher 
temperatures. 

The nonlinear thermal conductance is not sensitive to 
the Kondo correlations due to small contribution of low 
energy excitations to it. As a result the nonlinear thermal 
conductance n shows similar scaling behavior with S as 
the linear one. In other words, the nonlinear n can be 
calculated, with good accuracy, from its linear values. 

It is the nonlinear differential conductance which 
changes with S in qualitative way. In symmetrically cou- 
pled QD (5 — 1) it shows a maximum at zero bias, while 
the asymmetry tends to move the maximum to non-zero 
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voltages. All the other quantities change with S only 
quantitatively. 

V. CONCLUSIONS 

We have studied electrical and thermal transport in a 
quantum dot asymmetrically coupled to external leads. 
Those studies revealed that the most suitable tool to 
study the asymmetry in the couplings S is nonlinear dif- 
ferential conductance. The presence of the maximum at 
non-zero bias voltage directly signals that the system is 
asymmetrically coupled to the electrodes. In the linear 



regime electrical and thermal conductances show simple 
scaling behavior with S while thermopower does not de- 
pend on it. Finally, we also checked the Mott formula 
of thermopower and found that this quantity is a good 
approximation even in the presence of strong Coulomb 
interactions, provided the changes in the spectrum are 
confined to small region around chemical potential. 
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